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C/3 

O Abstract: We study the noise of the intensity variance and of the intensity 

'-(_> correlation and structure functions measured in Hght scattering from a 

5^ random medium in the case when these quantities are obtained by averaging 

over a finite number A^ of pixels of a digital camera. We show that the noise 

scales as 1/A^ in all cases and that it is sensitive to correlations of signals 

corresponding to adjacent pixels as well as to the effective time averaging 

(due to the finite sampling time) and spatial averaging (due to the finite 

H ' pixel size). Our results provide a guide to estimation of noise level in such 

Qh applications as the multi-speckle dynamic light scattering, time-resolved 

'""' correlation spectroscopy, speckle visibility spectroscopy, laser speckle 

,__! imaging etc. 

> 

r~^> OCIS codes: (030.6140) Coherence and statistical optics: Speckle; (110.6150) Imaging sys- 

QQ terns: Speckle imaging; (170.3880) Medical optics and biotechnology: Medical and biological 

^v/l imaging; (100.4550) Image Processing: Con'elators. 
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1. Introduction 

The statistical properties of optical speckle patterns resulting from scattering of light in a ran- 
dom medium are largely independent of the nature of the latter [ 1 1. Rearrangements of scatterers 
due to Brownian motion or flow lead to characteristic fluctuations of the speckle intensity with 
time. Probing these intensity fluctuations is a very common and highly sensitive non-invasive 
method to study the dynamics of complex fluids and biological systems. One of the earliest 
applications is dynamic light scattering (DLS), also known as photon correlation spectroscopy 
(PCS), where the fluctuations of the far-field speckle are analyzed 0. The method is widely 
used to study polymers in solution or for the sizing of sub-micron particles. The same method 
applied to turbid media is known as diffusing-wave spectroscopy (DWS). In both DLS and 
DWS, one usually makes use of the ergodicity of intensity fluctuations and replace ensemble 
averaging by time averaging. Both are typically applied to study fluctuations in the (sub-) mil- 
lisecond range. Thus, for a total measurement time of a few minutes, millions of fluctuations 
are sampled which provides an excellent signal to noise ratio. Using a single photon counter 
and a digital correlator, nanosecond time resolution is commonly achieved. Even though at the 
shortest times the signal to noise ratio is limited by photon shot noise |3|, in most practical 
cases this fundamental limitation is not the main concern and other sources of noise (such as 
sample purity, stability of the experimental setup, etc.) play the main role. 



Over the last decade, a wealth of new laser-speckle based experimental techniques have been 
introduced. Most of them are made possible by recent advances in optical sensor technology. 
A modem digital CCD (charge-coupled device) or CMOS (complementary metal oxide semi- 
conductor) sensor contains millions of detectors which allow massive parallel processing of a 
large number of signals corresponding to intensities of distinct speckle spots. The availability of 
an area detector has lead to numerous new applications in traditional far-field speckle detection 
such as multi-speckle dynamic light scattering El|5]|6l|2l[8l, time-resolved correlation spec- 
troscopy (TRC) |9|, speckle visibility spectroscopy (SVS) |10, 11] and has also enabled new 
speckle imaging techniques such as the near-field scattering (NFS) |12|, laser speckle imaging 
(LSI) 113, 14, 15, 161 or echo speckle imaging fTTlfTSl. 

The obvious advantages of massive parallel detection provided by an image sensor come 
at a price of specific and often critical limitations. The temporal resolution of digital cameras 
(typically, in the millisecond range) is inferior to the resolution of a traditional single photon 
counter such as a photo-multiplier (nanoseconds). In addition, to maximize the signal to noise 
ratio the experimental setup is usually designed in such a way that the size of an individual 
speckle spot is comparable to the size of the individual active area (pixel) of the camera. Hence, 
in contrast to the traditional PCS experiment, a camera pixel is not an ideal point-like detector 
and, furthermore, fluctuations detected by neighboring pixels are correlated. It would be, of 
course, desirable to increase the speckle size with respect to the size of the camera pixel (to 
mimic a point-like detector) and to reduce correlation of signals detected by neighboring pixels 
(to decrease noise in the statistical properties of speckle patterns determined by averaging over 
a group of pixels). However, these two objectives are mutually exclusive, at least if one wishes 
to exploit signals from all available pixels of the image sensor 

The limitations of digital cameras as image sensors appear most pronounced in such appli- 
cations as, for example, LSI. In LSI an image of a sample acquired by the camera is divided in 
a set of "meta-pixels", typically 25-100 pixels each. The spatially resolved information about 
statistical properties of the sample (an "image" in the statistical sense) is obtained by analyzing 
statistical properties of speckle (mean, variance, etc.) within each meta-pixel separately. Typi- 
cally, the measurements are made in backscattering from, for example, a biological tissue, and 
allow a full-field monitoring of dynamic properties, such as blood flow. The method is widely 
used for biomedical studies |fT9l l20l 1211 l22l |231 since it provides access to physiological pro- 
cesses in vivo with excellent temporal and spatial resolution. LSI is also becoming increasingly 
popular in soft material sciences as a probe of heterogenous dynamic properties 1 18, 26 1. Ob- 
viously, speckle imaging techniques require a tradeoff between spatial resolution and statistical 
accuracy |16|. Due to the rather low number of pixels forming a single meta-pixel (25-100 as 
already mentioned) the statistical accuracy and noise are of major concern for the data analysis. 

The goal of this article is to provide guidelines for a better understanding of the measured 
quantities and their fluctuations in speckle-based optical techniques for essentially all cases 
of practical interest. Whenever possible we provide analytical expressions that can be directly 
applied to the analysis of experimental data. In particular, we analyze the influence of a limited 
time resolution and provide an approximate treatment of spatial correlations between signals 
corresponding to neighboring pixels. We start by considering the variance of intensities in a 
stationary speckle patterns and then extend our analysis to the time correlation function and the 
intensity structure function of dynamic speckle patterns. 

2. Properties of stationary speckle patterns 

Consider coherent laser light scattered by a random sample in a typical light scattering ex- 
periment in reflection (see Fig. [Til. Usually the cross polarized channel is analyzed in order 
to suppress specular and low order scattering reflection. A digital camera provides us with A^ 
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Fig. 1. Experimental setup. A flat sample is illuminated with a polarized expanded laser 
beam. The diffuse reflected light is detected in the cross polarization channel by imaging 
the surface of the sample with a digital camera. The size of individual speckle spots in the 
image plane can be adjusted by changing the aperture of the camera objective. 



values of integrated (over the area of a pixel) intensity /«, a = l,...N, corresponding to A^ 
distinct pixels. A^ could be a (typically, large) total number of available pixels, but it could also 
be just a small subset of pixels carrying information about only a part of the optical field, scat- 
tered by a small part of the sample (see Fig. l2]i. Our purpose is to characterize the statistical 
properties of the optical field based on this information. Note that quite generally the statistical 
properties of the fully developed speckle patterns considered here are the same for measure- 
ments in the image plane or in the far-field fl^. Therefore our results should find applications 
in many research areas where speckle based techniques are employed. 

2.1. Negative exponential distribution of integrated intensities 

We start by considering a speckle pattern resulting from scattering of light in a solid sam- 
ple, where the positions of scattering centers are fixed. It is known that under quite general 
conditions, the intensity / of light scattered from such a random medium follows the negative 
exponential distribution: P{I)~{\/{I)) exp(— // (/) ) 1 1 1 . In the experiment we only have access 
to 



d\l{r). 



(1) 



pixel a 



which is the intensity integrated over an area a^ of a pixel of the camera. Here a — I,. ..,N 
indexes different pixels of the camera. As the first and the simplest example, we assume that la 
follows the negative exponential distribution as well and, in addition, that integrated intensities 
corresponding to different pixels are uncorrelated. The average value of intensity / and the 
variance of its fluctuations can be estimated from the data {/«} as 
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Fig. 2. Schematic representation of thie matrix of pixels of a digital camera. The full matrix 
is divided in a set of meta-pixels Nx x Ny = N pixels each. Spatially-varying statistical 
properties of the speckle pattern imaged by the camera (i.e., the variance of intensities c) 
are estimated by averaging over all pixels within the same meta-pixel. In our calculation 
(Sec. |2.3[ >, we are taking into account correlations between intensities at neighboring pixels 
in all directions. A given pixel (pixel 1 in the figure) have 8 neighbors: 4 neighbors of type 
2 and 4 neighbors of type 3. 



These statistical estimators are unbiased, i.e. the mean values of / and c found by averaging 
over an ensemble of realizations are equal to the actual values of the average intensity and of 
its variance, respectively: 



{{I-{I)f)^{I)' 



(4) 
(5) 



Here the angular brackets (. . .) denote ensemble averaging. The normalized variance {c)/{I)^ 
is equal to one in this case. The speckle contrast K can be defined as 



K 



c)/{I). 



(6) 



The values of / and c obtained in a series of measurements fluctuate around their means (|4| 
and (|5]). The variances of these fluctuations are 



or 



(c2) - (c)2 
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(8) 



where in the last line, the result after the "~" sign was obtained by taking the limit N ^t °°. 

To which extent are Eqs. (pUH]) relevant to realistic experiments? For a standard square matrix 
of adjacent pixels (see Fig. J2|l, two situations are possible depending on the relation between 
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Fig. 3. Left: False-color image of light intensities in a speckle pattern (60 x 60 pixels) 
obtained by using the smallest available aperture setting (//#=32). Intensity scale from 
to 10^ in arbitrary units. Right: Intensity correlation function obtained by fast Fourier 
transformation (FFT) of the raw speckle image (full frame 640 x 480 pixels). The data is 
quantitatively described by Eq. l|12[l with b = 1. 24a and fl = 1.09. 



the speckle size (i.e. the correlation range b of the scattered intensity) and the pixel size a. 
For small speckles (b <C a), different pixels are uncorrected, but the statistical distribution of 
integrated intensity la differs from the negative exponential one. In the opposite case of large 
speckles (b > a) the statistical distribution of /« approaches the negative exponential one, but 
intensities corresponding to neighboring pixels become correlated. The situation described in 
the present section — uncorrected pixels with negative exponential distribution of /« — can 
be reached by working with large speckles (b > a) but using only a subset of pixels (i.e., only 
pixels distanced by 1, 2 or even more pixels) in the analysis. However, using only a subset of 
all available pixels means that the total number of pixels is reduced. It is therefore desirable to 
extend Eqs. (pVISll to a more realistic model for the distribution function of la. 



2.2. Gamma distribution of integrated intensities 

For small speckles (b <C a) the integrated intensities /« coiTesponding to different pixels can be 
considered independent, whereas the distribution of each la, according to Goodman HI, can be 
modeled with the so-called gamma distribution : 



nia) = 



r(M)l(/> 



^a exp [ 



(9) 



Here the parameter /i > 1 depends on the pixel shape and size, as well as on the spatial corre- 
lation function g2(Ar) of intensity /(r): 
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where both integrations run over the area of the same pixel. If the sample is illuminated by a 
Gaussian beam, the intensity correlation function of the far-field speckle pattern is 



,^ ^ (/(r)/(r + Ar)) , 
^2(Ar) - ' ^ = 1 +exp 






(11) 




Fig. 4. Speckle parameter fl as a. function of the speckle size b divided by the size of the 
camera pixel a. Symbols: experimental results for the case of scattering from a solid sample 
(Teflon). Solid line: Eq. l |10| l. Dashed line: the approximate result in the small-speckle limit 
fc<Ca. 



whereas for a plane wave passed through a circular aperture. 



g2{Ar) = 1 



2Ji{Ar/b) 
(Ar/b) 



n2 



(12) 



In our experiments, the correlation function of intensity in the speckle pattern is controlled 



by the (circular) aperture of the camera objective, so that Eq. (12 1 is more appropriate for us. 
In both cases b quantifies the extent of spatial correlations of /(r) (i.e., the size of a single 
speckle spot). An analytic expression for jj. can only be obtained for very small speckles b<^a: 
jX ~ {a / b)'^ / {4-7t) . In Fig. Hwe compare the theoretical Eq. ( 10 1 with experimental results and 
find an excellent agreement. 

The expectations and the variances of / and c defined by Eqs. (J2]) and ([3]l can be calculated 
for the gamma distribution using the known expression for its statistical moments: 
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We obtain 
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where, again, the last result of the second line is obtained in the limit of A' ^- oo. As expected, 
these results coincide with Eqs. (pUHll for /i = 1, when Eq. (J9]) becomes identical to the negative 



exponential distribution. When jJ. > I, the variance of intensities c of the speckle pattern is 
reduced due to the effective spatial averaging over the area of a pixel that becomes comparable 
to the typical size of a speckle spot b [see Eq. ([T5|)]. At the same time, fluctuations of c from 
one meta-pixel to another are suppressed by a factor (1 + 3/jU)/4 < 1 [see Eq. (16 1]. 



Figure [3] (left) shows a speckle image of light scattered from solid Teflon taken for the case 
of large speckles that can be resolved in space (//# = 32). Be /(q) the Fourier transform of 
the image /(r). From the inverse Fourier transform of the power spectrum |/(q)P we directly 
obtain the non-normalized intensity correlation function gii^r) |[T1 [see Fig. [3] (right)]. In our 
example, for speckles of finite size b/a — 1.24, the intercept of the correlation function is given 
by l + l/jU = 1.92. 

2.3. Role of correlations between neighboring pixels 

The results obtained in the previous subsection can be extended to include correlations between 
neighboring pixels. In a matrix of adjacent square pixels (Fig. |2]) of side a larger than the corre- 
lation range b of the speckle field, it appears reasonable to take into account only correlations 
between integrated intensities la corresponding to neighboring pixels and to neglect correla- 
tions between intensities of more distant pixels. Note that each pixel has 8 neighboring pixels 
of two different types (pixels 2 and 3 inFig.|2]i. The impact of correlations on the average values 
and variances of / and c is described by two additional parameters 
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where, in contrast to Eq. (10 1, the two integrations in each equation run over two different 



pixels: the neighboring pixels 1 and 2 situated side by side in Eq. (17i or the pixels 1 and 3 



situated next to each other on a diagonal of a square matrix of pixels in Eq. ( 18 i, see Fig.l2] 

With this definitions in hand, we can readily calculate the average value of /, (/) = (/), and 
its variance: 
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The average value of c can also be readily found: 
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(19) 



(20) 



For A^ — > oo, Eq. (J20]) reduces to Eq. ([TSj. 

To obtain these results we gave special treatment to pixels situated at the boundaries or in the 
corners of the square matrix. These pixels have less neighbors, which is important when Eqs. 

I and ( [20| are used for small numbers of pixels A^ ^ 10. An interesting consequence of Eq. 

I is that the expectation value of c becomes A^-dependent. Because we saw in the previous 
sections that when different pixels are considered independent, (c) does not depend on A^ [cf. 
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Eqs. (Is]) and ( 15 1], we conclude that the dependence of (c) on A^ is a signature of correlation 
between signals corresponding to different pixels. The estimator (|3]l acquires a bias. Interest- 
ingly, the spatial range of correlations b or, more precisely, the ratio b/a can be estimated by 
observing the A^-dependence of (c) . 
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Fig. 5. Average variance of the intensity fluctuations (c) and its noise a^ as functions of the 
number of pixels N. Speckle pattern is recorded in the image plane for light reflected from 
a solid piece of Teflon. Exposure time is 1 ms, transport mean free path in Teflon /* ~ 0.25 
mm, camera pixel size a = 9.9 /im, magnification one. Left panel: analysis using a subset 
of pixels (all neighboring pixels omitted, full symbols) leads to a constant value of (c) [Eq. 
fT5|, dotted lines]. Analysis using all pixels (open symbols) is compared to the prediction 
of Eq. pOl (solid lines). Right panel: dotted and solid lines show predictions of Eqs. ^ 
and p.l\, respectively. The inset shows a comparison of the experimental values (symbols) 
and theoretical predictions foi H(fl) =Na^/{c) . Eq. 1 16 > (dotted line) andEq. ([2ll (solid 
line). 



The variance of c, C7^, can be calculated along the same lines as (c) . The calculation, however, 
is much more involved and rapidly leads to cumbersome equations. Instead of dealing with 
this lengthy and, anyway, approximate analytic calculation, we estimate a^ by a numerical 
simulation. We use a computer to generate a set of 4096 x 4096 speckle patterns with correlation 
lengths b — 2,A, and 8, respectively. These speckle patterns are then used to generate integrated 
speckle patterns with pixel sizes a varying from 2 to 64. The integrated speckle patterns model 
outputs of a typical CCD or CMOS camera. Finally, we numerically compute the dependencies 
of (c) and a^ on the number A^ of pixels, for fixed ratios b/a. For b/a <^ 1 the parameter 



H characterizing integrated speckle patterns [Eq. ( 10 1] is a function of b/a only. Hence, by 
inverting this function we can express b/a and hence a^/{c)^ as a function of ji instead of b/a. 
We know that Eq. (16 1 is valid in the limit of /i — ?► oo (or, equivalently, b/a -^ 0) which suggests 
that a^/{c)^ can be expressed as a series in 1/jti of which the terms of order {l/jJ.)'^ and (l//x)' 
are already known from Eq. ( [T6| . In the limit of large A^, the results of our simulations can be 
fit by adding a term quadratic in 1 //i : 
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where we introduced i/(/x) ~ 2(1 +3//! + 12//!^). This equation appears to describe our results 
quite well for < l//i < 0.8 (see Sec. |2.5| . For uncorrected pixels H{p.) ~ 2(1 +3/^) from 
Eq. ([T6]l. 



2.4. Properties of time-integrated speckle patterns 

If the scattering medium is not stationary like, e.g., a liquid suspension of dielectric particles, 
the intensity / of scattered light fluctuates in time. In such an experiment the intensity / of 
scattered light changes not only as a function of position r, but also as a function of time t: 
I = I{v,t). Fluctuations of / in time can be characterized by its autocorrelation function, 

giiT) - ^p , (22) 

where we omitted the position r from the arguments of /, for clarity. The analysis presented in 
the previous section applies to this situation too, provided that the sampling time T (i.e. the time 
during which the signal is accumulated to determine /«) is much shorter than the characteristic 
correlation time Tc of variations of/. In the opposite case, i.e. when T > Tc, the definition of /« 
should include not only the spatial integration over the pixel area, but the temporal integration 
over the sampling time as well. 

Under the assumption of small pixel size a -^ b and statistically independent /« following 
the negative exponential distribution, the effect of finite sampling time can be described by a 
parameter V lISJIfTTl : 

^ ^ ^ [82{^)-l]il-T/T)dx. (23) 



V T Jo 

The integrated intensity /« follows the gamma distribution (J9]l with v substituted for /z (l\. The 
average value (c) and its normalized variance are given by Eqs. ^5\ and ( 16 1, respectively. 



For pixel size a that becomes comparable to the correlation length b of the speckle pattern, we 
can take into account correlations between integrated intensities corresponding to neighboring 



pixels following the approach of Sec. 2.3 The outcome of our analysis is quite simple: all 



results of Sec. 2.3 — and, in particular, Eqs. (19 1, (20 1 and (21 1 — still apply but with jx. 



jX2 and /X3 replaced by /x x v, jU2 x V and /i3 x v, respectively. Here jx, jX2 and /i3 describe the 
impact of spatial correlations, whereas v accounts for time integration. The effects of spatial and 
time correlations thus fully decouple. This is due to the fact that the spatio-temporal correlation 
function §2(^1", t) — 1 = (/(r,f)/(r + Ar,f + t))/(/)^ — 1 decouples, in its turn, into a product 
of position- and time-dependent parts. 

2.5. Comparison with experiment 

In our experiment, linearly polarized light from a solid state laser (Verdi V5 from Coherent, 
wavelength A = 532 nm) is expanded and collimated to a beam of several centimeters in waist 
to create an approximately homogeneous illumination spot on the sample surface (see Fig.fT]). 
Liquid samples are contained in large glass containers whereas solid samples are measured 
in air The diffuse reflected light is monitored in the image plane in the crossed polarization 
channel with a CCD camera PCO Pixelfly (640x480 pixels of 9.9x9.9 \ivc? size, 12 bit) at 
magnification one. The estimated depth-of-focus of the imaging system is close to 0.1 mm. 
The actual size of speckle spots on the CCD chip can be adjusted by varying the aperture 
of the camera objective. For the simple configuration of a single lens placed at a distance / 
from the screen and a circular aperture with radius q the speckle correlation length is 1241 
b = 2l/q=[Al/kf]y.f/#. 

Our camera objective has an equivalent focal length of / = 50 mm and is placed at a distance 
/ ~ 100 mm from the CCD sensor. The /-number of the objective can be varied from //# = 2.8 
to 32. The accessible values of b are thus comparable or smaller than the pixel size b < a. 
The camera exposure time is adjustable with the shortest time used of 0.1 ms. The maximum 
speed of data acquisition at full resolution is 50 frames per second. The recorded images are 



corrected for slight spatial variations of the average intensity and, subsequently, the statistical 
properties are calculated from the recorded intensity values. We have carefully checked that this 
correction procedure does not alter the experimental results. Throughout this article we compare 
two different types of experiments with theory. In an idealized experiment we analyze only a 
subset of pixels where neighboring pixels are omitted. The procedure effectively eliminates the 
effect of spatial correlations between the pixels for all situations considered here. The second 
set of data is obtained by analyzing all available pixels, a situation typically encountered in 
actual applications. 

We first analyze the static speckle pattern of light reflected from a solid piece of Teflon 
(thickness «1 cm) with a typical photon transport mean free path I* ~ 0.25 mm. In Fig. l5] 
we compare our experimental results for (c) (left panel) and its normalized variance a^/{cy 
(right panel) with the predictions of the theoretical model developed above. The parameter /i 

is determined from a best fit to the data via the relation jx = (I) / (c) for N -^ °°. The values 
obtained for jx are found in excellent agreement with the theoretical predictions as shown in 
Fig. Here the parameter b for all settings has been obtained from the relation b °^ //# using 
the known value of b = 1.24 for //# = 32 (Fig.pb. For consistency, in the following discussion 
for dynamic media, we use the fitted values of jx for each setting. 

For the static samples, good agreement with theory is found for the A^— dependence of the 
mean and the variance of (c) . In particular, we clearly see that (c) is indeed independent of A^ for 
a subset of uncorrected pixels, whereas it acquires an A^-dependence for correlated pixels. The 
expected 1/A^ dependence of the variance of c is observed for both uncorrelated and correlated 



pixels, though with different prefactors. As follows from Eq. (21 1, the prefactor //(/x) depends 
on the degree of correlation between neighboring pixels that we quantify by the parameter jj.. 
It is shown in the inset of Fig. [3] The values of H following from the experiment are close to 
those expected theoretically. Small but visible discrepancies between data and theory in Fig.lS] 
are most probably due to the approximate nature of our theoretical model: even for uncorrelated 
pixels the statistical distribution of intensity does not follow the Gamma distribution exactly |[T] 
and weak correlations certainly exist not only between neighboring but between distant pixels 
as well. 

3. Noise in photon correlation spectroscopy 

3.1. Fluctuations of the intensity correlation function 

We now turn our attention to imaging and spectroscopy with speckle correlations. Instead of 
integrating the intensity of light scattered from a turbid sample with mobile scattering centers 



over a time exceeding the correlation time Xc of intensity fluctuations (as we did in Sec. 2.4 1, 
detailed information on temporal fluctuations of intensity can be obtained by choosing T ^ T^ 
and correlating speckle images separated by a time lag T larger than T . When intensities /«(? ) 
are measured for A^ pixels, the degree of intensity correlation can be estimated as 



N 

L 

a=l 



^(^) = ^ I [^« (^) - '(^)] [^« {t + ^)-i{t + ^)]. (24) 



where 



1 ^ 



'(0 = 77L/a(0 (25) 



N 



a=l 



is an unbiased estimator of average intensity (/); it is equivalent to / of Sec. [2] However, it will 
be important for us to keep track of the difference between /(f) and i{t + t) in the definition 
((24|ofc(T). 



The average and the variance of c(t) for independent pixels with la {t) distributed according 
to the negative exponential distribution are: 

(C(T)) = {lf[g2{T)~l], (26) 

<) _ g2(T)^(3-|)-2[g2(T)-l] ^ 1 ^ g2(T)[3g2(T)-2] 



(cW)' (A^-l)[g2(T)-lf A^ fe(T)-lf 

where the last expression is obtained in the limit of N ^ °°. 

As we saw in the previous sections, it is important to account for deviations of the intensity 
distribution from the negative exponential one as well as for correlations between neighboring 
pixels in the image, if a comparison with experiment is attempted. Because the intensity cor- 
relation function g2{Ar,z) — 1 = (/(r,f)/(r + Ar,f + t))/(/)^ — 1 decouples into a product of 
position- and time-dependent parts, we readily obtain 



(c(T)) = (/: 



1 4 4{VN-1) 



li ^2VNiVN+l) ^3NiVN+l) 



[g2(T)-l] (28) 



for the average of c(t). Here ji, fij and jj.^, account for spatial correlations of the instantaneous 
speckle pattern and are given by Eqs. ( 10 1, ( 17 1 and ([T8|), respectively. 



The calculation of the variance of fluctuations of c(t) with account for correlations between 
neighboring pixels appears to be much more involved. However, an approximate result can still 



be obtained by an ad hoc combination of Eqs. ( 2 1 1 and ( 27 1: 



<) „, //(M)^^ g2(T)[3g2(T)-2] ^ ^29) 



(c(t))2 8A^ fc(T)-l] 



To compare Eqs. (26l-(29i with measurements we have carried out a time correlation ex- 
periment on a slowly relaxing sample. We dispersed about 3.3% of polystyrene microspheres 
(diameter 710 nm) in an aqueous solution of cetylpyridinium chloride/sodium salicylate (100 
mM CPyCl/60 mM NaSal). The resulting surfactant solution strongly scatters light with the 
transport mean free path /* ~ 73 /zm at 532 nm. A detailed characterization of the system is 
given in Ref. Il27l . It displays strongly viscoelastic properties with a slow terminal relaxation. 
For our measurement we keep the sample at room temperature T ~ 22° C which sets the re- 
laxation time to several tens of seconds. We first determine the fall autocorrelation function 
from a combined photon PCS and CCD-camera based experiment |7|. This yields the tempo- 
ral autocorrelation function g2(T) of scattered light (see the inset of Fig. l6]l. We then perform 
measurements of the noise of speckle correlation coefficient c(t) for different sizes of the cam- 
era objective aperture corresponding to the values of jx ranging from 1.29 to 5.15, and as a 
function of the number of pixels A^. In all cases we observe that the normalized variance of 
c(t) scales with 1/A^ as expected from the theory. In Fig.p^we show N(J^/{c)^ as a function 
of the time lag T. The general rule is that the normalized variance of the correlation coefficient 
c(t) increases with T. Experiment and theory are in fairly good agreement if the analysis is 
performed on a subset of pixels (right panel of Fig. l6]), especially at small jj., when the statis- 
tics of intensity fluctuations is close to negative exponential. The agreement is worse when all 
pixels are analyzed (left panel of Fig.l6]l, even though the theory describes the general trend of 
our data quite well. This might either suggest that we were either unable to eliminate residual 
experimental errors or that additional factors, not taken into account in our theoretical model, 
enter into play when important correlations between signals corresponding to different pixels 
are present. 
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Fig. 6. Noise of the speckle correlation coefficient as a function of the time lag T (symbols). 
Solid line shows prediction by Eqs. {19\ (left) and \27\ (right). Inset of the left panel shows 
the intensity correlation function g2('^)- 



3.2. Fluctuations of the intensity structure function 

Instead of the intensity correlation function g2{T), fluctuations of intensity can be characterized 
by a structure function 



D{T)^m)-i{t+T)Y 

where for the case of Gaussian statistics 



cmt). 



.„ = ffl£W(i±l)ft=.M«)-«M 



(30) 



(31) 



The intensity structure function (ISF) is a more direct measure of the dynamic activity and has 
several advantages as compared to the intensity correlation function (ICF) ||3][T8!|. While both 
quantities are directly related in the limit of perfect measurement statistics, the ISF is known 
to outperform the ICF in accuracy when the collection time is limited. Furthermore, the ISF is 
less sensitive to low frequency noise or drifts |3|. An unbiased estimator of D(f) is 



1 '^ 
*(T)-T. Y, [Ia{t)-Ia{t + ^)f. 
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For the negative exponential distribution of integrated intensities Ia{t), assuming no correla- 
tion between intensities at different pixels, we obtain 



(..(T))2 
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(33) 
(34) 



Accounting for finite pixel size a that can become comparable to the correlations length b of 
the speckle pattern, we obtain 
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(35) 
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Fig. 7. Noise of the of the speckle structure coefficient as a function of time lag T (symbols). 
As predicted by the theory, the noise is independent of T. Lines show predictions of Eqs. 
l |36[ l (left) and ( |34^ (right). The inset of the right panel shows the intensity structure function 



Once again, the calculation of the variance of fluctuations of ^(t) taking into account correla- 
tions between neighboring pixels appears to be much more involved. However, an approximate 



result can still be obtained by a (though ad-hoc) combination of Eqs. (21 1 and (34i; 
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A comparison between Eqs. (34i, (|36]) and experimental data is presented in Fig. l7] The 
theory correctly predicts that the noise of .s(t) scales with 1/A^ and that it is independent of T. 
However, theory and experiment do not show quantitative agreement: theoretical lines lie either 
higher (blue and red lines in the left panel of Fig. 17]) or lower (all other lines in Fig. 17]) than the 
data. The reasons behind this disagreement are the same as in the case of intensity correlation 
coefficient c(t) in Sec. 3.1 



4. Conclusions 

In the present paper we presented a study of noise in modern speckle correlation and imaging 
techniques. The noise originates from replacing the ensemble averaging (assumed in theoret- 
ical models) by averaging over a finite number A' of pixels of a digital camera. In particular, 
we studied the noise of the speckle intensity variance c, and of the intensity correlation and 
structure coefficients c(t) and ^(t), respectively. The variances of all these quantities decrease 
as 1/A^ when the number of pixels is increased and depend in complex ways on the spatial (due 
to the finite pixel size) and temporal (due to the finite sampling time) correlations of intensity 
in the speckle pattern. For stationary speckle patterns, we obtained quantitative agreement be- 
tween measurements and the theoretical model that we developed in this paper. For dynamic 
speckle patterns, theoretical predictions reproduce general trends of our data but fail to provide 
a fully quantitative description. We believe that this is due to the approximate character of our 
theoretical model (neglecting correlations between distant pixels, gamma distribution of inten- 
sity at a single pixel) as well as to imperfections of our experiment (impossibility to achieve a 



perfectly uniform illumination of the sample, necessity to work with weak signals, etc.). De- 
spite the absence of quantitative agreement between theory and experiment in the latter case, 
the results presented in the paper provide an important starting point for estimation of the noise 
level in such applications as the multi-speckle dynamic light scattering, time-resolved correla- 
tion spectroscopy, speckle visibility spectroscopy, near-field scattering, laser speckle imaging 
and echo speckle imaging. 
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